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it Introduetion 


In a nuelear power reaetor neutrons emitted from the fission proeess of uranium 
255 must be slowed down inside of the reactor by a moderator if they are to propagate 
the fission reaetion. Beeause of te importance in reactor design of knowing the 
spatial distribution of these slowed neutrons many studies have been made of the 
slowing down proeess. Probably the earliest suesessful method of treating the slow- 
ing down theory ean be attributed to Amaldi and Fermi. In their model, eollisions 
between neutrons and the nuelei of the moderator are treated as free partiele eol- 
lisions . 

ne has ealeulated that the effeet of ehemieal binding of the moderator 
atoms ean be negleeted in the treatment of the eollisions between neutrons and modera- 
tor atoms. The Fermi slowing down theory should therefore be valid for the slowing 
down of neutrons whose energy is less than the binding energy of a moderator atom in 
its lattiee. The purpose of this experiment is to test the validity of the Fermi 
theory when applied to neutrons whose energy is less than the ehemical binding of 
the moderator. 

Graphite as a moderator is ideal for this purpose beeause the earbon atom is 
bound to the graphite lattiee with an energy of about 20 ev. In all other respeets 
graphite conforms well to the limitations imposed by the slowing down theory. 

A monoenergetie souree of neutrons whose energy is between 20 ev and a lower 


limit of 1 ev imposed by restrietions to be explained later, was not obtainable . 


A A <3 


», the spatial 


A monoenergetie souree of neutrons of about 24 kev was used, therefore 
distribution of 4.8 = l ev (4) gold resonanee neutrons was obtained and eompared by 
the Fermi theory to the 1.44 ev indium resonanee neutron spatial distribution obtained 


by Roberts, Hill and MeCammon $®? 











EE Theory 


It may be of interest to review the assumptions and limitations imposed in the 
derivation of the slowing down theory. Consider a fast neutron that makes an elastis 
eollision with a stationary unbound nueleus and transfers a portion of its energy 
to the nueleus. The energy lost by the neutron depends upon its initial energy, the 
angle through whieh it is defleeted and the mass of the nueleus. If the target atom 
is heavy, earbon for example, and the kinetie energy assosiated with it seonsists of 
tke thermal energy, its veloeity ean be neglected when eompared to the veloeity of 
a lev neutron. “ Further, if in a center of mass system the probability of the 
neutron being scattered in any small solid angle A, n is a constant, it may be show 
that the average loss of energy by the neutron depends only on its initial energy and 
the mass of the target nucleus. The average loss of logarithm of the energy per col- 


lision, y » is independent of the initial neutron energy and is defined as and—is— 


equal to 





id ~ bn, - Lnk, = 1+ kom (a) 


Ae on c 
where M = Am) 
and 


EB) = initial energy of the neutron 


Gy 
DO 
! 


energy of the neutron after collision 


> 
i 


mass of the target nucleus. 


When fast neutrons are released at a constant rate from a point in an infinite 
medium in which the above described collisions occur, the neutrons will diffuse away 
from the source and slow down until they are in thermal equilibrium with the medium. 


If the sicwing down is treated as a continuous process and the number of neutrons lost 














by absorption during slowing dowmm is negligible, the number of neutrons that slow 


dowm past energy E in one second fer om® ata distance r from the source is 


given by 


re 


R% = 
q(r.E) = [ert (ey) >” 6 4°C (8) (2) 


Qo = source strength number of neutrons per sec. 


q(rs,E) = slowing dowm density i.e. mumber of neutrons slowing down past 
energy E per sec cm 


E) 
1 


2 
E ¢ 8 (1 SA 


é 


T(E) = Fermi age = a (3) 


=. = NG. = macroscopic cross section for scattering 


N = number of target nuclei per om® 
GO. = microscopic cross section for scattering 


BE, = initial neutron energy 


E, = energy of neutron at detection. 


iby ome is independent of neutron energys T (BE) may be integrated to give 


T (8) ” —— on « KBr ra (4) 
32. F (1-5) d d 














Treating the slowing down of neutrons as a continuous process rather than a dis- 
continuous process introduces the largest error when one attempts to find the slowing 
down density for age zero or at distances far from the source for any energy above 
thermal cee e227) 

Because the slowing down density is related to the neutron flux, foil activation 


measurements can be used to determine the slowing down density. The relation is 


given below: 


B) 
Q(r,E)dE = _a(rB) (&, 
=,F5 
where Q(r,E) denotes the neutron flux and is identically equal 
to n(r,E)v 
n(r,E) = number of neutrons per em? in energy range d& 


v = scaler value of the neutron velocity cm per sec. 


For the purposes of this experiment, foil activations determined at various 
distances from the source were plotted, and [°(£) was determined from the slope. 
The 4.8 * 1 ev neutron age can be compared to the age of 147 em* for 1.44 ev neutrons 


determined by Roberts, Hill and McCammon by equation (4) 


Bo Au 
Cau ~ Cig 7 Ae BD 


Om - 147 = 15.425 Un +t (6) 














III Equipment and Techniques 


Except for the gold foils and counting instruments the equipment was the same 


as used by Roberts, Hill and McCammon and much of the following data were obtained from 


their report 


lo 


Zo 


(5) 


Graphite colwm. 

The 6' x 6' x 11' colum shom in Fig. (I) was composed of 4" x 4" 
graphite stringers, the central stringers along the vertical axis were 
removable for inserting the foils. Its size was such that it was effective- 
ly infinite in size for the neutron ages investigated. The average density 
of the graphite was 1.652 g/cm®. The scattering cross section of the carbon 


temas! 10K ne 


Neutron source. 

Practically monoenergetic neutrons were obtained from a source con- 
sisting of about 1 curie of antimony inside of a spherical shell of beryl- 
lium. The shell had an overall diameter of about 2.5 cm and weighed 9.24 
grams o 

Because the antimony has a rather high capture cross section for neutrons, 
a depression of the flux might be expected within a couple of mean free paths 
from the source. The absorptions of its owmm neutrons may also affect the 
half life of the antimony. This change is shom in the appendix to be less 
than 0.0067% and will therefore not affect the 60 day half life of the source 


appreciably. 














UNCLASSIFIED 
DWG. 9147 
FOIL HOLDER 


SOURCE 





FIG. I 
GRAPHITE COLUMN 














5. Detecting Foils. 
Gold foils measuring 4 x 6.55 cm and 103.1 mg/om* thick were used to 
detect the 4.8 ev neutron flux. The maximum difference in the weight was 
6%. During the exposure, the foils were enclosed in cadmium boxes shown in 


Fig.(II) to shield them from thermal flux. 





Two thin wall glass Geiger Mueller tubes enclosed im a 2" thick lead chamber 

to reduce the background counts were connected through preamplifiers to scales of 64. 
The tubes were normalized with activated gold foils over a range of counting rates 
that varied between 5,000 to 15,000 counts per minute. The normalization factor was 
found to be constant over this range of counting rates. 

Frequent counts were taken of two long half life foil standards to ascertain 
that the counting efficiencies of the instruments were constant. The two standards 
used were a Bay ae emitter with a half life of one year and a natural uranium 


‘Genta 








UNCLASSIFIED 


OWG. 9148 
FIG. 1 
CADMIUM BOX 
FIG. 


FOIL HOLDER CYLINDER 

















A 


Neutron Flux Measurements 

The isotopic neutron flux was detected in the following manner. Foils protected 
from the thermal flux by cadmium boxes were inserted at a number of positions along 
the principal axis of the graphite colum. During the exposure time, which consisted 
of about three days, the foils were prientediin the plane perpendicular to the 
principal axis of the colum. After this activation time they were removed from 
the colum and from the cadmium boxes, wrapped inside of aluminum cylinders, shown 
in Figo III and position on the Geiger Mueller tubes by a jig shown in Fig. IV, 
and counted for about eight hours. Immediately after one side of the foils had 
been counted, they were inverted and a similar count taken from the other sidee 
By use of the normalization equation given in appendix 3 the counts obtained from 
the foils were converted to activity based on a 21.65 day exposure. The calculated 
activity of both sides of the foil was averaged to eliminate the current component 
of the flux. It is this averaged activity that is proportional to the isotropic 
flux(5s7) o These data are given in coluwm 7 of Table I. 

To evaluate the condition of the graphite colum, indium foils were similarly 
exposed at a few positions in the columm and countede The calculated activities of 
the foils were compared with the activities reported by Roberts et al for these 
positions. <A tabulation of these measurements are given with their data in Table 


II. 
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TABLE II 


Measured saturated activity of indium 
Distance of foil from foil covered with cadmium 
neutron source ; 


Check Points| Roberts »Hilly McCammon 


room + -gam | Act(x,2) x 1.0912) Act(r,E) |LogioAct(r,E) 


3.551 12.61 3590 3.54720 
6.091 | 37-10 3236 3.50001 
8.631 Th.4Q | 2924 3.45271 
11.171 124.79 3.40870 
geal 187.96 3.35017 
16.25 264.06 3.28007 
18.79 353.06 3.19846 
21.33 45h .97 3.10584 
23.02 529.9 3-05030 
23.82 567 4 1030 3.02168 
28.10 789.6 2.85205 
28.90 835.20 2.81713 
33.18 1100.9 2.61363 
33.98 1154.6 365.0 2.57262 
38.26 1463.8 2.33646 
29 0% 1524.1 2.295h2 
43.34 1878.4 2.05423 
48.42 23k 5 1.71675 
49.27 2427 25 1.67822 
53.50 2862. 3 1.42570 
1. 33245 














Experimental Results 

A plot of the logarithm of the foil activity versus radii squared is 
given in Fig. V. A Gaussian distribution, equation 2, in these coordinates 
will appear as a straight line with a slope of - 1/4T( E). 

The neutron age may be obtained quite reasonably as 128 ft 2 cm from the 
slope of the central portion of the curve. Near the origin much of the devi- 
ation from equation 2 may be attributed to the higher energy absorption 
resonances in gol ? © At a radius of about 22 cm this effect becomes 
negligible and the curve conforms quite well to a simple Gaussian out to a 
radius of about 45 cn. 

The experimental difference between the indium resonance and gold 
resonance neutron ages is 19 + 2 cm. By the use of equation 6 the calculated 
age difference between 4.8 + .1 ev neutrons and 1.44 ev neutrons is 18.6 = .3 
om“. These results then indicate that within the Beerary of this experiment 
the Fermi slowing down model is valid for neutrons slowing down in graphite 


from 4.8 ev to 1.44 ev. 


VL Corrections. 


The use of a spherical neutron source and finite size detector foils instead 
of a point source and point detectors ae. effect om the results of this ex- 
periment. Roberts, et al ; have shown that the effect of the spherical neutron 
source is to change the constant of proportionality between the foil activation 
aud the neutron slowing dow density. They also stated that a similar analysis, 
which is given in Appendix 2, would demonstrate that the only effect of the finite 


foils is to again change the proportionality constant. Because the absolute value 


of the neutron slowing dow density is not of importance in this experiment, the 
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corrected constant was not evaluated. 

Unfortunately foils become activated by neutrons whose energies are not with- 
in the resonance emergy range, i.e. resonamce energy minus to resomance energy plus 
ome half the average emergy loss per collision. To correct for this non resonance 
activation it is required that the variation of capture cross section with neutron 
energy be known. 

Such a correction was attempted following the analysis given im Appendix 4 
and using for the capture cross section for gold the Breit Wigner formula 


: ES Goi 
Sg) a EB [@ - £) * 2/4] 


where E, = 4.8 ev, the neutron energy for resonance peak 


Fa 


0481 ev, the nuclear energy level width 


Oo = #700 barms, the capture cross sectiom at resonance peak. 


This procedure resulted im a corrected age for the gold resomance 


neutrons of 127 cm“. Because the values of [7 is perhaps much smaller than 


0481 ev, the corrected age should perhaps be somewhere between 127 cm’ and 128 cm*. 
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TIL APPENDIX 


le To show that the absorption of neutrons affects the antimony half life a 
negligible amount, the assumption is made that every neutron produced by the 
source is self absorbede Wattenberg \®) states that approximately one half of 
the antimony disintegrations produce the 1.8 Mev //raye In the same article the 
=28 


cross section of the y (Be-n) reaction is given as 9e7 x 10 om”. 


The approximate rate of forming neutrons is given by 


NO pT 
ry (1-9 Be abe Be, 


® thickness of Be shell 





A s decay constant for antimony 


= 
Pw 
i 


@ |rate of decay of the antimony | 


= 4317 x 107°* Be atoms per cm® 


oy 
i 


Or # 0.4565 cm, thickness of Be 
aBe 


N «().000134. 
rate = A (a2 ) = 0.000067 NX. 


The actual rate of decay of the antimony may be obtained from 





aN = 
s- =~ ANG S. (1 = 9702000134, 


” 


dN 
—— =~ AN(1 = 0.000067) 
dt 
and 
x apparent # (] ~ 02000067) A 
2e To show that no foil geometry corrections are necessary, the following assumptions 
are madée 
&e Point neutron source, proof that the actual source appears as a point source 


is given in ORNL~201. 














be Circular, thin foilse 


Ge The spatial resonance flux distribution is Gaussian. 


radius of foil 


Me) 
ag 





T @ thickmess of foil 


r 8 distance of foil from source 


=x2/V (E) 
¢ (87) =a = __%° 


§ ye Ey , é a B, [at x (BQ)? 
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-r°/4% (B,) 


4U(E,)17 Q. No; T _ . R 
| dA = oe ° | af AX (Bo) r 2/4 U(Eo)a P 
b = 


vol ED e,(are (g,))9/? 


4% (E,)7 Q, Nox T 


EQ z,far v (e,)]°/ 





Act = B E = eRe (Bo) | ere /4t (Ey) 














Since R® and t (E,) are to be constant the total foil activation varies with 
. : ar” /4.% (BE) 
distance from the source as oe °% 
Se All foils were corrected to the activation expected if they had been placed in 
the colum the first day, exposed for 21.9 days and counted immediately. 


N, # number per om? of Aut97 


N # number per cm® of active gold atoms 


A, = antimony decay constant & 0.0155/daye 

r»> @ gold decay constant ® 0.2567/daye 
The source strength varies with time as 

y(t) =a "ies 
If the number of Aut98 atoms formed is small compared to the number of Avl97 
originally present,the rate of formation of Aul98 may be expressed as 


~ =s- A, NtNo, $(t) 


but Ne To is considered constant ;therefore let A NeTa = B 
= A,% 


GN = 
+A, N=Be 


The boundary condition is that N = 0 when ¢t # 0, 
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The number of Au nuclei formed at any time t, is then 


ny =——B___ ~ Arty - 
i = ~ v 
Ne 7 Ay 
When the foils are removed from the colum,the gold decays as 


dN 
et -NA, 








ee <=: “eee 





The number present No at the end of this cooling interval t, will be 


-Azt, «= Agt, Ag te 
é e de 


The counts obtained, C, in the counting interval te will depend upon the counter 
efficiency A 
tz 


- Agt 
cs-d #2. 21,2, | a at 


0 


© 


- Ax =Aot, Ast ~Ao% 
le 1 By Re 2°, 


The maximum number of counts that could be obtained would occur if t 1 were equal to 


_&n Ao/ >, 
7 


given for antimony and gold, t, for maximum activation in the colum is 2109 dayse 


» tg equal to zero and t, equal to infinity. For the decay constants 


C = a_B e = @ 


~21.9), “21.9, 
max ae ) 
We = Ay = ~ 


By forming the ratio of maximum counts to the actual counts and multiplying by the 
gold decay constant the activity normalization factor for 21.9 days exposure time 


and zero cooling time is obtainede 


=F 0 ) ~21.9 7 
nex * “Aye, T AQT Agt Agts 
(e - 6 


1s "aes ) 


- 











To correct the activity to that which would have been attained if the foils had 


been inserted into the colwm on the lst day, Act. ax must be multiplied by the 


= Xt 
source factor 1/e aa 


oe The tq is the time elapsed between the insertion of the 
first foil on the first day and the insertion into the colum of the foil one 
wishes to normalizee 


The normalizing equation becomes 


= “21.9 A, a0) | 
2 © # © 


lst day) (e” Aye, | af dott) Ast (a:eeelh Azts) wv Arta 


oad Te 


4- To account for low energy gold activation the following analysis is appliede 


d/l = differential solid angle 
u = cos @ 
ds = differential area 


ds 
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Assumptions : 


Be 


be 


Ceo 


de 





Isotropic fluxe 


The scattering cross sections of the shield and foil materials can be 


neglectede 


Absorption of neutrons of energy E does not perturb the neutron flux at 


any lower energye 


Plane geometryo 


Consider the foil at position ro 


(1) 


(2) 


(3) 


(4) 


Y (E,r,)dE = number of neutrons with energies between E and 

E + dE which cross unit area per seconde 

Y (E,r,)dE u = # number of neutrons with energies between E 
and E + dE which cross unit area normal to the foil per second 





in the solid angle d+ 3 277 du. 
“MN, 7742 u 
o (E,r,)dB u i< e = the number of these neutrons which 





penetrate the shield to a depth of z cm per second if Ty is the 
absorption cross section and N, is the number of absorbing atoms /cm® 
of shielding materiale 


Of the neutrons that penetrate this far the number absorbed in 


layer dz om thick per second is 





= on a/a 
(Er )dButt , ** Neds 
47¢ u 
or as 277 du Sav 
-N,O7,z/a 


Y(B, r,)daBuSee + * No-dz 
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AOO= sfu 

(5) ¢(B, r,)B u -| 7% (-—2 dz) 
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2 
absorbed in thickness Tj in angle due 


-N, 0 ,T./u 
(6) Y(B,r,)u = dee } iy/ 


between E and E + dE getting through one om’ of shield per second 


“N07, T./u 
2 9 (Byr,) we fa eo i Bec 


= number of neutrons with energies 





with angle du about ue 
(7) By applying similar reasoning to the neutron absorption in the foil 
and multiplying the result by a factor of 2 to account for activation 


from the other side the activation becomes 


N_o-_T /u Nto-'t'/y 
Y(Esro)u e 1 iil lee —_ dE du 


where N'q-'T' refers to the gold foile 
(8) Since Y (E,r,) is independent of u, the equation in (7) may be 
integrated over u with limits of 0 to le By a change of variable 


f= J AL e@ @ 
wu = 
-—~- the activation becomes 


20,000 ev 20,000 ev 
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£ 2 Zz . i ee 
Cd Cd 
cut off cut off =— 
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where ~ E, — = —— 
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OF and o-' are functions of energy 


o/s (8) 


[en ecny] § 5 2 


> and E to be evaluated for graphitee 
Ss 
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p (Bsr) = 


(9) The correcting factor used is 
_ (468 0634) ev 
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(408 = 0.34) ov 
( 20,000 ev 


Cd cut off 
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